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We present the first computation of the thermodynamic 
properties of the complex su(3) Toda theory. This is possi- 
ble thanks to a new string hypothesis, which involves bound 
states that are non self-conjugate solutions of the Bethe equa- 
tions. Our method provides equivalently the solution of the 
sm(3) generalization of the XXZ chain. In the repulsive 
regime, we confirm that the scattering theory proposed over 
the past few years - made only of solitons with non diago- 
nal S matrices - is complete. But we show that unitarity 
does not follow, contrary to early claims, eigenvalues of the 
monodromy matrix not being pure phases. In the attractive 
regime, we find that the proposed minimal solution of the 
bootstrap equations is actually far from being complete. We 
discuss some simple values of the couplings, where, instead of 
the few conjectured breathers, a very complex structure (in- 
volving Ea, or two Eg) of bound states is necessary to close 
the bootstrap. 

PACS: 72.10.-d, 73.40. Gk 



Integrable quantum field theories (QFTs) based on 
su(2) have been the subject of intensive studies over the 
past many years. New theoretical tools like the quantum 
Q-operators Q, the Destri De Vega and Klumper and 
Pearce equations H, the connections with spectral de- 
terminant theory |3j, or the relations with elliptic curves 
and duality j| have revealed mathematical structures of 
remarkable depth; they have also made possible the com- 
putation of quantities of experimental interest - and typ- 
ical of strong interactions - in various low dimensional 
condensed matter systems 

Formal developments, as well as practical applications, 
would largely benefit from an extension of these results 
to the case of other Lie algebras, in particular su(n). 
The situation here is somewhat embarrassing, however. 
Although the pillars of the su(2) case - the XXZ chain 
and the associated sine-Gordon model - have been under 
control for a long time, even the simplest sit(3) case is 
only very partially understood. 

One of the difficulties here - and, from the field theory 
point of view, one of the most interesting issues at stake - 
has to do with unitarity. Indeed, the simplest integrable 
generalizations of the sine-Gordon model are the complex 
[] affine su(n) Toda theories defined by the Lagrangian: 
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where A > 0, at, ■ ■ ■ , form the root system of the 

classical lie algebra a n _i, a Q = — X)j=i a i * s ^ e ne £ a ~ 
tive of the longest root. The conformal weights of the 
perturbing operator in ([l]) are A = A = In the 

following, we shall parameterize A = ^-f-- The theory 
described by (0) is obviously non unitary at the classi- 
cal level. The fascinating possibility was raised HQ that 
it could nevertheless describe a unitary field theory in a 
sufficiently strong quantum regime. This possibility was 
ruled out in the interesting paper , and we confirm and 
extend their observations here. 

From a practical point of view, unitarity is not such 
a key issue. In fact, the most interesting applications 
of complex Toda theories are potentially found in disor- 
dered systems of statistical mechanics, where Toda the- 
ories based on superalgebras naturally seem to appear 
leading most likely to even stronger violations of 
unitarity. More crucial then are the questions of com- 
pleteness of the bootstrap, the physical meaning of the 
bound states, and the calculation of physical quantities. 

The main progress in the study of complex Toda theo- 
ries have been based on non perturbative S matrix anal- 
ysis, following the pioneering work of ||. One of the dif- 
ficulties in this approach for su(n) is the appearance of 
a large number of poles in the S matrix elements, whose 
signification is not entirely clear: it was argued, after 
careful analysis of several cases, that most of these poles 
were not physical, and occurred rather by mechanisms 
generalizing Coleman Thun's [[TT|| . 

The issue of the completeness of the bootstrap in |||7]] 
could be settled by a study of the thermodynamics, and 
a computation of the central charge in the UV, using the 
thermodynamic Bethe ansatz (TBA) p^] ^. However, 
for the imaginary affine Toda theories (or, equivalently, 
the corresponding anisotropic quantum spin chains), the 
TBA has never been written so far, because of the com- 
plexity of the set of solutions: no natural "string hypoth- 
esis" had been proposed, up to now. 

In this letter, we present the first solution of this vex- 
ing problem in the case of su(3). This allows us, in par- 



1 Real Toda theories involve entirely different issues. For a 
recent review see Iptl. 



2 The approach developed in Q bypasses some difficulties in 
the study of the thermodynamics, but does not give much 
information on the spectrum itself [J13|. 
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ticular, to show that the quantum theory denned in (|l]) 
for n — 3 is never unitary, even in the strong quantum 
regime, and that it presents considerably more bound 
states than expected. 

Our main technical progress is an understanding of the 
solutions of the Bethe equations for systems based on 
su(3). Except in the exactly symmetric case, these so- 
lutions (the equivalent of the string hypothesis for the 
XXZ chain) had never been found. We hope that their 
understanding will spur new development in the area: 
generalizations of the works mentioned in the introduc- 
tion, as well a calculations of physical properties in the 
super Toda case, seem particularly timely. 

The problem we want to tackle is largely equivalent 
to the su(3) generalization of the XXZ spin chain. This 
integrable chain has been known for a long time |l4[ ,and 
reads 
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Relations of this system with the quantum field theory 
are two-fold. At the physical level, there is a simple in- 
tegrable perturbation - obtained, in the quantum inverse 
scattering framework, by introducing heterogeneities of 
the spectral parameter, which amount to a staggered in- 
teraction - giving rise to ([!]) in the continuum limit, as 
can easily been shown using bosonization, or other ar- 
guments ]l5| ]. In that correspondence, 7 = j. At a 
more formal level, observe first that the scattering ma- 
trices proposed in |l6[ ] are non diagonal. To study the 
thermodynamics of the gas of excitations in the S matrix 
approach, one needs to write wave functions, and im- 
pose their periodicity. This condition involves "passing" 
a particle through a set of other particles with which it 
scatters non diagonally: the phases obtained in this way 
are, in the inverse quantum scattering framework, eigen- 
values of a monodromy matrix []l7f . This monodromy 
matrix is nothing but the transfer matrix associated with 
the hamiltonian (||) (actually, it involves a slight general- 
ization of this hamiltonian with a mixture of the funda- 
mental representation and its conjugate), with however a 
renormalization of the anisotropy parameter, which be- 
comes 7 = The case 7 = ^, for instance (the equiv- 
alent of the XX chain, which is not a free problem here 
however), is a lattice regularization of the Toda theory 
for A = i; it also turns out in the monodromy problem 
for the Toda theory at A = |. 

This double correspondence is very useful. If one is 
able to solve say the field theory problem for a particular 



value of A 



-, one should also be able to diagonalize 



of (@) for 7 = j in turn means knowing the spectrum of 
the monodromy matrix for this value of 7. This is nothing 
but being able to solve the Toda theory for A = -^-l 

Rather than dwell into the lengthy technical details, 
we would like to discuss first the solution of the lattice 
model for 7 = |. We consider the slightly more gen- 
eral hamiltonian, where there is an arbitrary mixture of 
the fundamental representation and its conjugate. The 
lattice Bethe equations are 
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the lattice model for 7 = j , since it is just a discretization 
of this quantum field theory. But knowing the spectrum 



where we have not written a crucial but complicated 
sign on the right hand sides, and the energy reads E = 
— 2V — 2V — -, — . A combination of numerical 

t—i cosn z cosh y 

studies and analytical arguments led to the identifica- 
tion of the following sets of roots in the thermodynamic 
limit. The y and z can both be real, or both have an 
imaginary part equal to 7r (antistring). In addition, it is 
possible to have complexes, with a two string z centered 
on an antistring y, that is: y = r + in, z — r ± i-| (here r 
is a real number, and we identified r + 3z-| and r — z-|), 
and the same thing with y and z reversed. The existence 
of these complexes is easy to understand. Suppose for 
instance that z has a positive imaginary part, so the lhs 
of the second equation in (|J) blows up as N3 — > 00. It is 
then necessary to have the rhs also blow up, which can 
be accomplished if there exists a y such that y = z + i?. 
The same goes if z has a negative imaginary part, so the 
complexes proposed are the minimal possible structures 
leading to real Bethe equations. Together with the real 
and antistring solutions, they reproduce the 3+3 degen- 
eracy expected for the fundamental solitons. In addition 
however, another type of complex (which we call yz in the 
following) is possible, of the form z = r-H^f , y — r+i^f, 
together with the conjugate. Usually, one would expect 
that these two complexes actually come glued together 
to ensure reality of the Bethe equations and the eigen- 
values, resulting in a sort of "quartet" (l9). This does 
not seem to be the case here. Rather, to reproduce the 
correct entropy or central charge, one needs to treat their 
densities as independent. 

These complexes are the ones which dominate the ther- 
modynamics. As usual, there are many other solutions to 
the Bethe equations. The existence of solutions which are 
not invariant under complex conjugation being rather un- 
usual, we illustrate it briefly. Figure 1 shows the yz com- 
plex obtained by a numerical solution of (3) for iVg = 
and for different values of N3. In this example - which 
seems the simplest possible - the half sum of the imag- 
inary parts of y and z goes to 7r in the thermodynamic 
limit, but their difference does not go to §. We conjec- 
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ture that this would however be the case for the over- 
whelming majority of such complexes: we checked that 
only in this case, are the correct scattering theory and 
thermodynamics recovered. 



12 si 
15 si 

13 si 
21 si 
24 si 
27 si 
30 si 
33 si 
36 si 
39 si 



FIG. 1. A numerical solution of (3) containing a yz com- 
plex for iVg = for different values of Nz = 6, 9, • ■ ■ , 39. For 
39 sites we show all roots belonging to this particular state, 
whereas for all other lattice lengths only the yz complex is 
shown. 

The lattice model without heterogeneities has a con- 
tinuum limit which is a conformal field theory, made of 
two bosons compactified on a triangular lattice. The lat- 
ter is unitary, and we have checked numerically that the 
imaginary parts of the eigenenergies of (^) scale to zero 
faster than their real parts, ensuring reality of the con- 
formal weights, and, actually, unitarity of the conformal 
field theory, despite the non-hermicity of the hamilto- 
nian. Meanwhile, the lattice model with heterogeneities 
does have a continuum limit described by a complex Toda 
theory, with A = |. In this case, the complexes corre- 
spond to physical particles, the bound states discussed in 
HQ . The fact that they are described by solutions of the 
Bethe ansatz which are not self conjugate corresponds to 
the shocking fact that, in the quantum field theory, their 
S matrix is not a pure phase (and "naive" unitarity, in 
the sense of the scattering theory, is broken). 

Accordingly, a crucial feature related with the yz com- 
plexes is that they give rise to non real kernels in the 
continuum Bethe ansatz. For instance, the kernel for the 
scattering of a real z through a yz complex z — r + i ^f , 

y _ „ i „• 5tt ;„ 1 d - 



if is 



ln$ = -l/cosh(z-r + if). Al- 
though at equilibrium the densities of the two types of 



conjugate complexes are equal, it is necessary to let them 
vary independently to get the correct entropy. The phys- 
ical results will be the same as the ones that would be 
obtained with a different theory, where the scattering ker- 
nels would be real, and given by the real part of the true 
kernels. This is what we will consider in the following, 
to make the notations simpler. 

The lattice model whose equations are (||) occurs also 
in the problem of diagonalizing the monodromy matrix 
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The complexes just discussed have no phys- 



ic + if:, and the same 



ical meaning then: they just correspond to pseudo par- 
ticles of zero mass. However, the fact that they are not 
self conjugate means that eigenvalues of the monodromy 
matrix will not, in general, have modulus one. Contrary 
to early claims 0, unitarity is strongly violated in the 
theory, even in the absence of bound states. The physi- 
cal implications of having "monodromies" which are not 
pure phases are not clear to us: presumably S matrix 
elements have to be considered as formal objects used to 
build wave functions, and cannot be given a reasonably 
meaning in terms of scattering processes [J 

This basic structure generalizes easily to the case t in- 
teger. For 7 = j, the Bethe equations contain in addi- 
tion a term of interaction between y roots, and between 
z roots. The solutions are the usual y or z real, anti- 
strings, or 2, . . . , t — 1 strings, plus y-t strings centered 
on z antistrings, and the same thing with y and z re- 
versed. In addition, the yz complexes now are of the 
form z = r + iir — if:, y = r 
with y and z exchanged. 

To proceed, one has to use these results to solve the 
monodromy problem. At A = ^-f-, passing a soliton 
through a set of N^,N^ solitons at various rapidities 
gives rise to a"'phase" which is an eigenvalue of the mon- 
odromy matrix. The latter is essentially the transfer ma- 
trix associated with (||), but for 7 = jf^. Its eigenvalues 
are obtained by generalizing slightly the Bethe ansatz 
equations, so the arguments in the lhs of (||) are shifted 
by the corresponding rapidities, and using the general- 
ized string hypothesis described previously. 

A few technical details are involved, which are com- 
pletely equivalent to what happens in the su(2) case [ ^p| . 
We will only describe the end result, which is quite sim- 
ple. For the case t — > 00, the TBA has been known for 
a long time [^l|j22[ , and has a structure that mimics the 
known one for su(2), with an infinity of massless nodes 
corresponding to the usual strings of the Bethe equations. 
Like for the su(2) case, the introduction of anisotropy 
truncates this to a finite number of strings: the trunca- 
tion has to be completed by the proper "end structure" 
of the diagram. Like in the su(2) case, this structure, 



A somewhat related problem occurs for instance in the 
massless description of conformal field theories. 
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for su(3), is given by the Bethe roots for 7 = g. After 
some manipulations, the results are the following. We 
call a, a h the densities of the solitons at rapidity 9, and 
m the mass of the fundamental soliton. Then 



a 3 + a 3h = m cosh 9 + * (pf + cr 3 ' 1 ) 



(4) 



and similarly for 3, with pq = <r. We recognize here 
the standard equations for the minimal model |l6|| . In 
addition, in this untruncated case, we need the "closure" 
relations, which are the key to the whole problem. They 
read here 

p?_3 + Pt-3 = <f>* (pt-4 + Pt-2 + pl-3 + Pi, + Pl 2 ) 

+ i>*(pb+Pb) 
Pt-2 + Pt-2 =<t>* (Pt-3 + Pt-2 ~ Pi, - Pi 

-1> *(/% + /%) 



pI+pII 



Pt-2 + Pt-2 



(5) 



and similarly for 3; finally 

pi +p 3 b h = V*(p?-3 + P?-2 



p 3 ai ~ Pl 2 



<P*(pI + pI) 



(6) 



In these equations, the kernels are defined by their fourier 
transforms 4> — 2 eo 1 sh - and ip = c ° s ^\!x ■ subscripts 
di, i = 1,2 stand for antistrings or t strings centered on 
antistrings; b stands for yz complexes. 

We set e- £ °/ T = a/a h , e~ e "/ T = p^/pi, n = l,...,t- 
2, e~ <L "/ T — pa,i/Pa i e~~ tb / T = Pb/Pb (color labels are 
kept implicit here). Introducing the usual variables x a — 
e -e a /T ^ -j- ne TBA in the UV limit reduces to the system 
(x t _ 2 = l/x a ): 

a+*i) i/2 (i+^r /2 

X„ = (1 + X„_ 1 ) 1 /2(l + Xn+1 )^ (l + 

X t . 3 = (1 + X t _ 4 ) 1/2 (1 + X a ) 3 / 2 (l + ^) _1/2 (1 + X b ) 



-1/2 



(l +2 ; t _3) 1/2 (1 + ^ 



-1/2 



(l+X a )- 1 (l+^)~ 1 



x b = (l + x t -3) (i + ^J (1 + aj o r 2 (H- Si)" 1 
The solution of this system is 



We define />) = i J f(6)e- mx9 ^ dO. 



x b 



_ (j + l)(j+4) • _ n 
— 2 ; J — u ) - 

t+1 



,i-3 



4t 



(8) 



In the IR limit, we get an identical TBA, but with the 
replacement t — ► t — 1 exactly as in the SU{2) case [ pp[ , 
because the first two nodes become infinitely massive, 
60 = 00. The whole TBA system is actually quite simi- 
lar to the one for the sine-Gordon model: the "left part" 
can be conveniently encoded in a ladder diagram with 
base the Dynkin diagram of ci2, and in the isotropic limit 
t — * 00, this is all that matters. For finite t however, the 
diagram has to be closed to the right, and the closing 
terms are more complicated than for sine-Gordon: in ad- 
dition to 3+3 nodes standing algebraically for the repre- 
sentations 3 and 3 (in the sine-Gordon case, there are two 
such nodes corresponding to the self-conjugate funda- 
mental representation), the closure requires the nodes as- 
sociated with the non real solutions of the Bethe ansatz, 
whose algebraic meaning has eluded us up to now. The 
TBA diagrams for SU(2) and SU(3) are shown in Fig. 2. 
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FIG. 2. The TBA diagrams for ST/ (2) and SU(3) 
Using the basic identity p3]: L(x) + L(y) = L(xy) 



L 



x(i-v) 

l-xy 



L 



y(i-z) 

l—xy 



, one finds 



(7) 



t-3 
3=0 

+2L 



3 + 5.7 + 4 



P + 5j + 6 



6L 



t - 1 



2t 



t + 1 



= (*-!) 



(9) 



and therefore, c = 2(t— 1) — 2(t— 2) = 2, as required. Note 
that for the sine-Gordon case, the sum of dilogarithms for 
the UV (IR) diagram was (t — l) 2 ^- (resp. t — 2) instead, 
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generalizing nicely the combinatorics to the case of rank 
two. 

Still like in the sine-Gordon case, the same TBA dia- 
gram supports the scattering theory for generalized su- 
persymmetric extensions: by putting the mass term on 
the k th node, the central charge now becomes c = 
in agreement with the central charge for SU (3)& theories. 

The Toda theory can be twisted by coupling a charge 
to the total diagonal U(l) current. The kinks come in 
the form of 3 and 3, and the total charge can be recast 
in terms of the solutions of the associated, non diagonal 
scattering problem, as Q oc (t — 1) J /of_ 2 — + (3 — > 3). 
The thermodynamics is more complicated to carry out 
than in the sine-Gordon case however, densities of holes 
and particles mixing up in a complex fashion in the right 
hand side of the universal Bethe equations. It can still be 
shown that, with a fugacity —1 for the solitons, the dia- 
gram truncates, leaving only the nodes j = 0, . . . , t — 5, 
for t > 5, the well known result for the minimal model 
SU{3)i x SU(3) t -4:/SU{3)t-3- Note that more nodes dis- 
appear than in the sine-Gordon case: the case where the 
minimal theory has c = now corresponds to a dimen- 
sion A = | in the Toda theory. 

The study of the points A = ^=i, t integer, therefore 
presents few surprises, except maybe the non unitarity, 
due to the complex solutions of the Bethe equations. We 
do not expect a qualitatively different behaviour when t 
is not an integer, provided A lies in the repulsive regime, 
A > 3/5, where no bound states are expected. In con- 
trast, the attractive regime turns out to be considerably 
more complicated. 

From a technical point of view, the points A = j , could 
be studied as the continuum limit of the chain with 7=7 
still, but with the opposite sign of the hamiltonian. This 
observation is crucial in the case of the XXZ chain, and 
allows one to essentially use the same type of solutions of 
the Bethe equations to study both the attractive and the 
repulsive regime: it is well known that the TBAs for these 
two cases are very closely related. In the sm(3) case, the 
situation is more involved: depending on the sign of the 
hamiltonian, it is not the same type of solutions of the 
Bethe equations that determines the thermodynamics. 
This can be illustrated in the case 7 = §• Taking say a 
chain made only of 3's, the ground state is y and z real for 
one sign, y antistring and z real for the other. This still 
is like for the XX chain. However, the other solutions of 
the Bethe equations are different: for instance, instead of 
having two string z centered on antistrings y's, one has 
now two string z centered on real y's. Things get worse 
as t increases, and solutions that were possible but did 
not contribute to the thermodynamics for one sign of the 
coupling, start being important for the other. 

The net result is that we have not been able to repro- 
duce the simple c = 2 result for any point in the attractive 
regime except the A = i case. To illustrate the nature 
of the difficulties, let us concentrate for a while on the 



theories for which A = — 1 is of the form A = 2 " 3 +1 , 
or A = g~ . In that case, the deformation parameter 
for the quantum group symmetry of the scattering the- 
ory q — — exp(i7rA) is a cubic root of unity, and the 
q-dimension of the 3 or the 3 vanish exactly, (3) q = 0. 
This means that the RSOS truncation HJfH of tne the_ 
ory is very simple: the kink part entirely disappears from 
the scattering, and one is left with a diagonal scatter- 
ing for the scalar bound states. This scattering leads to 
straightforward TBA calculations, which one can com- 
pare against the predictions of conformal field theory. 
The disagreement with the scattering theory proposed in 
M is considerable. Let us give two examples. For A = |, 
the analysis of predicts only a pair of bound states: 
to reproduce the effective central charge c e g = 6/7 to- 
gether with its apparent conformal weight A„g = 1/7, 
one needs in fact 6 bound states, essentially reproducing 
the well known E§ scattering theory |^(|. Note that in 
this theory, there are pairs of conjugate particles, but also 
self-conjugate ones: the issue of whether bound states in 
su(3) should always appear in pairs is not clear, and this 
counter example demonstrates it is not always the case. 

If the appearance of Eq is somewhat surprising, an even 
bigger surprise is encountered for A = |. In this case, 
the effective central charge is c e g- = 1, and the apparent 
weight A e g = 1/16. The scattering theory in predicts 
two pairs of bound states, i. e. 4 scalar particles. To 
close the bootstrap, it turns out that one needs in fact 16 
particles, with a scattering theory that is a sort of double 
Eg, p7j : if one labels the Eg particles by j, the su(3) 
particles by Bj, then one has the identity S B .§ j S BkBk — 
Sjk- this guarantees that the central charge in the UV is 
2 x ~ = lQ, as required. 

The ratio of the actual number of bound states and 
the minimal one is actually growing very quickly, and 
it is not clear what to do to close the bootstrap even 
for the simplest non trivial reflectionless point, A = |, 
where the conjectured scattering theory ]^,|28[|, involving 
the fundamental solitons A ,A , excited solitons A\,A\, 
and breathers Bx,B\, £?2,i?2, gives only c = 1.90821. 
Note that this point is quite special, since the particles 
can be gathered into 2 eightuplets of particles of mass 
M and y/3M, indicating the presence of a hidden 
symmetry ]29] , |3"o| ] . Introducing more particles increases 
c, but we haven't found a closure of the bootstrap, even 
after introducing large numbers of bound states []. 

The repulsive regime for su(n), n > 4, presents similar 



5 This doubling is similar to what happens in the scattering 
theory for the 3 state Potts model, with c = 4/5, compared 
with the scattering theory for the Yang Lee model, with c„g- = 
2/5. 

6 Let us stress here that the possibility of having missed some 
CDD factors in the original scatterings should be considered 
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difficulties: the scattering theory of @^8| is incomplete 
even at the point A = i, and there again, more bound 
states have to be introduced. 

In conclusion, except for the simplest case of sit(3) in 
the repulsive regime, the scattering theory of the com- 
plex Toda theories is not under control, in our opinion. 
The large body of literature on the subject (see e. g. ]3^ ] 
for some recent developments) has identified a minimal 
structure for it, but reality is more complicated (a fact 
well known to some experts). How complicated is not 
clear to us: without an algebraic understanding of the 
bound states, tackling the general case is a very labori- 
ous task, and it is not even clear that the bootstrap will 
close with a finite number of bound states, except for 
some special values: clearly, more work is needed in this 
direction. 

Meanwhile, the solution of the su(3) theory in the 
repulsive regime opens the way to several problems of 
physical interest, in particular study of su(3) anisotropic 
Kondo problems, and su(3) tunneling problems, on which 
we hope to report soon. 
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